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“ FOREWORD 


Program WIND (Winds In Non-Uniform Domains) was a cooperative 
applied research program conducted by the USDA Forest Service 
Reoye anc, the U.S. Army, and supported by other public and private 
PeOoperacors. The purpose of Program WIND was to obtain data for 
the evaluation of computer-based predictive mesoscale 
meteorological and particulate dispersion models. The program 
was conducted in northern California over an area approximately 
80 km x 80 km bounded on the south by Sacramento, on the north by 
Redding, on the east by the Sierra Nevada Range foothills, and on 
Macewest by Interstate 5, and consisted of four phases = Summer 
Phase 1, Winter Phase II, Spring Phase III and Fall Phase IV. 

The field work began in the summer 1985 and ended in the fall 
e9e7.., The USDA Forest Service participated by conducting 
esctralcewake and dispersion studies during Phases I, III and IV. 
This was the first meteorlogical mesoscale study of this scale in 
North America in terms of the type, number, and frequency of 
meteorological data collection. Cooperation among the 
participating agencies and the advancement in meteorological 
instrumentation made this study possible. 


The FS participated actively in Program WIND by providing a 
program director (Jack Barry, Washington Office, Forest Health 


ProeecetOnmotatt), a public affairs officer (Ann Westling, Tahoe 
NattOnal Forest), and assistant program director (Bob Ekblad, 
Missoula Technology Development Center). The FS provided 


administrative management, coordination with the public and news 
media, and contracting and leasing of real estate for the 
numerous study and sampling sites. 


The program approach was to characterize atmospheric conditions 
by obtaining meteorological data during each of valley wind 
patterns associated with the four seasons and during different 
times of the day. Data were used to evaluate and enhance a 
hierarchy of predictive meteorological models. Prior to Program 
WIND the models had suffered for need of a database extensive 
PuouommlOecharaccerize and quantify the physical processes at 
play in the atmosphere. These include models and integrated 
Mogelewthac predict (1) wind flow over a mountain-valley complex, 
(eminGluenceror forest and crop canopies on wind flow, (3) 

Ot eers10n Ofesmoke, air pollution; dusts, seeds, pesticides and 
other agricultural materials and (4) effects of aircraft wake on 
Of pens toneCre agricultural and forestry sprays and dry materials. 


CoeelT cally thes os interest focused on models which predict 
dispersion of smoke from controlled burns and movement of air 
Pollutants, spray drift, penetration, and deposition of sprays 
into and on forest/range canopies. Data collected during Program 
WIND was also used by the FS to develop and to make FS spray 
dispersion models simpler and easier to use in the field. 





The U.S. Army was primarily interested in a capability to predict 
winds as they affect obscuration from smoke and dust over an area 
menging from one km to 6400 km. 


While the FS reported most of its findings from Program WIND at 
symposia, in professional journals, and in FS reports, some data 
escaped reporting. This report and others being issued under the 
Forest Health Technology Enterprise Team (FHTET) and Missoula 
Technology & Development Center (MTDC) report series covers some 
of those unreported studies conducted by the FS in conjunction 
with Program WIND. By providing reports within the FHTET and the 
MTDC report series, the data will be more accessable to a broader 
audience through FS and other library searches. 


Key Words: 


forest meteorology, spray models, dispersion models, aerial 
application, forest canopy, spray drift, smoke behavior, forest 
pest management, program WIND. 


Honn W. Barry Harold) Thrstle 
Director, Forest Health Program Leader 
Technology Enterprise Team - Davis MTDC 
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EXECUTIVE SUMMARY 


On October 15 and 26 twenty-eight passes were made by F-15 and F-16 high speed fighter 
aircraft normal to a row of ten towers holding propeller anemometers. These anemometers 
measured the wind speed generated by the aircraft wake at 32 locations. The vortex circulation 
time history was then inferred by a generalized vortex-locating algorithm developed for previous 
comparisons with slower speed aircraft data from Program WIND. These data results were then fit 
by least squares to the AGDISP/FSCBG circulation decay model. Representative atmospheric 
turbulence levels were recovered for these runs, and the decay coefficient was determined. 


These results extend the applicability of the Program WIND results to the higher aircraft 
speed regime. The quantity defined by decay coefficient times turbulence level is shown to be 
larger than the previous Program WIND results, and suggests that a mean value of 0.74 m/sec is 


an acceptable input parameter for higher aircraft speed regimes. 
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1. INTRODUCTION 

Program WIND ("Winds in Nonuniform Domains") was performed in phases jointly by 
the United States Department of Agriculture Forest Service and the United States Army at several 
sites in northern California. Phase I was conducted from May through July, 1985 at a cleared and 
forested site (Foresthill Seed Orchard) and an almond orchard (Hennigan Almond Orchard north of 
Chico); and Phase III, from late April to early May, 1986 at a cleared sloping site in the Sierra 
Nevada foothills near Red Bluff. In both of these experimental studies, anemometer tower grids 
recorded the ambient vertical velocity time histories as various slow speed aircraft repeatedly 
traversed normal to the grid. These digitized velocity traces produced an aircraft wake signature 
that could be used to infer the strength and lateral and vertical motion of the aircraft vortex pairs 
generating the traces. A previous report (Ref. 1) summarizes the examination of this data base by 
using a generalized algorithm to locate the trailed vortex pairs in all available data, and to infer their 
_ decay properties in the atmosphere. This decay effect was then quantified as a decay coefficient 
range for subsequent input in the AGDISP and FSCBG codes (Refs. 2 and 3). 


In the Phase I tests the tower grid included anemometers for both horizontal and vertical 
velocities, with a total of 32 anemometers available for the analog-to-digital system employed to 
record the data. This study determined that the horizontal velocity signals would be easily polluted 
by any crosswind velocity component, giving a decidedly incorrect impression of the vortex 
position and strength. Also, the closeness of the tower grid (the towers were spaced 2.7 to 6.1 
meters apart with seven towers) meant that no more than 20 to 40 seconds of data could be taken 
before the vortices would drift off the grid. Thus, in the Phase III tests the towers contained only 
vertical velocity anemometers, and were spaced uniformly apart at 6.1 meters (with ten towers). 
A longer sampling time was used to track the vortices across this wider tower grid. 


On October 15 and 26 a series of flights at Dugway, Utah, were made by high speed 
fighter aircraft over a tower grid similar to the Phase III tests and previous high speed data (Ref. 
4). The towers were extended to 15.2 meters in altitude and spaced 12.2 meters apart (again ten 
towers) in an attempt to capture the passage of the vortex pair. The sampling time was continued 


until the anemometers appeared visually to cease rotating. 


The test procedure followed at Dugway is summarized in Section 2 of this report. Section 
3 reviews the generalized algorithm used to reduce the data, and Section 4 compares the vortex 


circulation decay data with previous results. Conclusions are offered in Section 5. 


2. SUMMARY OF TESTS CONDUCTED 
The test procedure on-site at Dugway, Utah, is reviewed in this section of the report. 


251 Test Aircraft 


The aircraft used in all passes over the tower grid were either an F-15 or F-16 high speed 
fighter aircraft. The intent of the test series was to extend the aircraft flight speed well beyond the 
aircraft speeds in Program WIND and confirm the results of previous tests (Ref. 4). 


Qe2 Test Site 


All tests were performed over the flat terrain near Dugway, Utah, with the tower grid set 
perpendicular to an access road (Tango Road) to offer a visual cue to the pilot. Initially, smoke 
pots were also used to identify the road. Ten towers were spaced uniformly apart and 
instrumented with propeller anemometers as shown in Figure 2-1. Data acquisition was housed in 
a command trailer centered on the anemometer grid approximately 10 meters north of the tower 


line. The site was the same as for the previous tests with an F-4. 


The towers used to mount the anemometers were telescoping masts extended to their 
maximum height of 15.2 m. The towers were placed in a single line normal to the road at 12.2 
m intervals, yielding a grid span of 110.0m . The anemometers were mounted vertically at 3.3, 
8.3 and 13.8 m from the ground on each tower. In addition to the vertical sensors, U-V-W 
sensors were mounted on tower 0 at 14.1m andontower9 at 8.3m. 


Dies Anemometers 


Gill anemometers manufactured by the R.M. Young Company, on loan to the U.S. Forest 
Service from the Transportation Systems Center of the U.S. Department of Transportation were 
used for this study. Four-bladed propellers (19 cm in diameter) were coupled with the DC 
generators to complete the anemometer and were calibrated at 1800 rpm. The anemometers were 
mounted on the towers oriented vertically to measure the vertical velocities of the wake. The 
anemometers were electrically connected to the data acquisition system with filtering capacitors in 
accordance with manufacturer's recommendations. The wake velocity signals were sampled and 
recorded by a digital data acquisition system during the aircraft flight over the towers. U-V-W 
sampling was also done prior to each aircraft run. 
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2.4 Data Acquisition System 


The data acquisition system consisted of an IBM Portable PC with two Data Translation 
DT2801 auxiliary boards to sample and digitize 32 channels of analog voltage signals from the 
anemometers. Sampling was carried out at a rate of 100 samples/sec (each anemometer was 
sampled every 0.34 sec). The full scale analog voltages +1.25 volts were converted to the 
digital representations 0 to 4095, and stored in memory. A maximum of four minutes of data 
could be sampled continuously (before reaching memory limits) and stored for post-test conversion 
to engineering units and further analysis. Sampling began before the aircraft flyover, to provide 


U-V-W data to extract the mean wind conditions. 
2.5 Test Matrix 


The objective to the test program was to collect high speed aircraft wake data (vortex 
trajectories) to compare with the slower speed results found from the Program WIND tests. The 
complete matrix is given in Table 2-1. No data was collected for runs 10 and 27 because data 
reduction was not complete before the aircraft returned for the next pass. Run number 24 was 
skipped. 
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NY NY NY WH WH HN HN HN HW NHN KR YR FR HH Re 
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DAY 


10/15 
10/15 
10/15 
10/15 
10/15 
10/15 
10/15 
10/15 
10/15 
10/15 
10/15 
10/15 
10/15 
10/15 
10/15 
10/26 
10/26 
10/26 
10/26 
10/26 
10/26 
10/26 
10/26 


10/26 
10/26 
10/26 
10/26 
10/26 


TABLE 2-1. DUGWAY TEST MATRIX 


TOWER AIRCRAFT 


Center 
Center 


East 
6-7 

Center 
Center 
Center 
Center 
Center 
Center 
Center 
Center 


Center 
Center 


Center 
Center 


F-16 
F-16 
F-16 
F-16 
F-16 
F-16 
F-16 
F-16 
F-16 


F-16 
F-16 
F-16 
F-16 
F-16 
F-15 
F-15 
F-15 
F-15 
F-15 
Ek 
F-15 
F-15 


F-15 
F-15 


F-15 
F-15 


2-4 


HEIGHT 
(Meters ) 


30.5 
30.5 
30.5 
30.5 
30.5 
30.5 
30.5 
30.5 
30.5 


30.5 
30.5 
BOs) 
30.5 
30.5 
OES 
30.5 
30.5 
30.5 
3055 
S10ks) 
30% 
30.5 


30.5 
B02) 


30.5 
30.5 


SPEED 
(Knots) 


300 
300 
300 
400 
400 
500 
500 
550 
550 


350 
350 
300 
300 
300 
300 
300 
300 
400 
400 
500 
500 
550 


450 
350 


450 
450 


REMARK 


Trial Run 
Empty Tanks 
Full Tanks 
Empty Tanks 
Full Tanks 
Empty Tanks 
Full Tanks 
Empty Tanks 
Full Tanks 
No Data 
Empty Tanks 
Full Tanks 


No Tanks 
Empty Tanks 
Trial Run 
Full Tanks 
Empty Tanks 
Full Tanks 
Empty Tanks 
Full Tanks 
Empty Tanks 
Full Tanks 


Number Skipped 


Empty Tanks 
Full Tanks 
No Data 

Full Tanks 
Empty Tanks 


3. GENERALIZED ALGORITHM FOR DETERMINING VORTEX TRAJECTORIES 


Reference 1 summarizes the generalized algorithm for inferring the location and strength of 
the vortex pair traveling through an anemometer tower grid. That analysis is reviewed here. 


An aircraft wake may be represented by a simple vortex pair (with its image pair below the 
surface, Figure 3-1), each vortex of which may be characterized by a velocity field of the form 


ep a (3-1) 


2n (Re ” R?) 
Cc 
where 


= velocity magnitude 
= vortex circulation strength 


radius from vortex center 
c = vortex core radius, a constant 


rep el tes| os 
II 


The resultant velocity at any point in the flow field (in particular, at an anemometer 
location) is then the algebraic sum of the velocity contributions from the four vortices acting on the 
flow. The aircraft vortex pair is located in a coordinate System relative to the tower grid 
(Figure 3-2). The wake model introduces four unknowns that must be deduced from the data in 
any run at any time at which data are collected. These unknowns are the vortex circulation strength 
and the three spatial dimensions 


h = the height of the vortex pair above the surface 
the semispan of the vortex pair (the half distance between their centers) 


the offset distance of the aircraft centerline from the leftmost anemometer 
tower 
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Aircraft wake physics are only approximated by the simple velocity law given in Eq. (3-1). 
Thus, the intent of the generalized analysis is to seek a solution for these four parameters in any 
Tun at any time so as to minimize the error in the velocity predictions at all of the anemometer 
locations. With four unknowns and many more anemometer locations, the problem may be cast 
into a least squares analysis by defining an error E as 


N 
Eee (Wn- w,)s (3-2) 


n=1 


where the overbar denotes the data, and the index n is from 1 to N , the total number of 
anemometers. The error E_ is a positive definite quantity. The vertical velocities w are 
determined by summing the contributions of the four vortices in the model, calculated for specific 


values of the parameters T , h , s and d . The data w are found directly from the test results. 

A crucial physical observation is made at this point: if the experimental data were 
normalized by any positive value, the resulting velocity vectors could still be used to infer the 
spatial position of the vortex pair in the anemometer field irrespective of the actual vortex 


circulation strength. This observation implies that locating the vortex pair involves the solution for 


h , s and d without regard to the value of TI . 


The approach used here normalizes the vertical velocities in the error equation by their 


respective root mean square velocities to give 


r- (33) Gy) : 


where 


fa 2 

Wrms = wee Wh (3-4) 
— 1 are 

Wrms — Vn Se il 





Because the predicted vertical velocities are all linear in vortex circulation strength I , this 
modification to the error removes I from the equation for E* and recasts the problem into one 


involving the three parameters h , s and d. 


Equation (3-3) for E* is nonlinear in h , s and d _, but may be analytically differentiated 
with respect to these unknowns to give 


F, = dE*/dh 
F, = dE*/ds (3-5) 
Fy = dE*/dd 


Clearly, when all three derivatives in Eq. (3-5) tend to zero, the slope of E* also tends to 
zero and E* reaches a local minimum. These F functions are nonlinearly related to h , s and 
d . At any value set of (h,s,d) , the local values of the F's in Eq. (3-5) may be determined 


analytically to generate partial derivatives of the form OF/dh , OF/@s , OF/dd to construct the 
model equation for F (representing Fh, Fy and Fy ) of the form 


Ey OE ie 2 BE is (3-6) 


oh Os od 


The partial derivatives in this equation are known as influence coefficients. For any value 
of (h,s,d) , the values of the three partial derivatives and F may be generated. Since the solution 
requires F to tend to zero, the left-hand side of Eq. (3-6) may be specified by 


Abeer (3-7) 


The complete system may be seen to be a three-equation system for the values (Ah,As,Ad) 
to add to the present values (h,s,d) to give new values (h+Ah,s+As,d+Ad) to recompute E* 


and drive the solution to a minimum. This iteration procedure is followed at every time in the test 
run being analyzed. 


| Once values for (h,s,d) have been determined, the original error E may be differentiated 
to give 


Fp = 0F/oF = 0 (3-8) 


to find the minimum of E and determine I . This equation is solved easily because the equation 
for F is linearin I . With the solution in hand for (1\h,s,d) , these values are used as initial 


guesses for the next time in the run being analyzed. All results have been computed with a core 
size of R, = 0.15 meters. 


Figures 3-3 to 3-28 present the results of applying the generalized algorithm to all 
applicable data. In all cases the time of 0 seconds is the time of initiation of the data acquisition 
system. The algorithm was started after 5 seconds of data collection in an attempt to begin data 
reduction when a vortex pair could be expected in or near the tower grid. The horizontal or Y 
distance is measured from the center of the tower grid 4-5. The vortex pair track horizontally and 
vertically is shown in these figures along with the inferred circulation strength. The dashed curve 
in the circulation plot is a best fit to the vortex decay model. It may be seen that most of the runs 
provide excellent correlation with the circulation decay model of Eq. (4-1) below. 
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Figure 3-4: Track for run 2 at Dugway. 
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Figure 3-5: Track for run 3 at Dugway. 
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Figure 3-6: Track for run 4 at Dugway. 
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Figure 3-8: Track for run 6 at Dugway. 
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Figure 3-7: Track for run 5 at Dugway. 
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Figure 3-9: Track for run 7 at Dugway. 
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Figure 3-10: Track for run 8 at Dugway. 
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Figure 3-11: Track for run 9 at Dugway. 
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Figure 3-12: Track for run 11 at Dugway. 
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Figure 3-14: Track for run 13 at Dugway. 
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Figure 3-15: Track for run 14 at Dugway. 
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Figure 3-16: Track for run 15 at Dugway. 
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Figure 3-17: Track for run 16 at Dugway. 
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Figure 3-18: Track for run 17 at Dugway. 


DISTANCE (CM) 


HEIGHT (MM) 


EStMAx274SEE) 


qQ2X 


1) Is) 30 4B 50 50 
i6Q_18 


12 
8@ 


46L 





0)2) 68 


28 38 40 
RIMS GSi203 


Figure 3-19: Track for run 18 at Dugway. 
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Figure 3-20: Track for run 19 at Dugway. 
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Figure 3-21: Track for run 20 at Dugway. 
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Figure 3-22: Track for run 21 at Dugway. 
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Figure 3-23: Track for run 22 at Dugway. 
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Figure 3-25: Track for run 25 at Dugway. 
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Figure 3-26: Track for run 26 at Dugway. 
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Figure 3-27: Track for run 28 at Dugway. 
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Figure 3-28: Track for run 29 at Dugway. 
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4. DETERMINATION OF DECAY COEFFICIENT 


The vortex strength decay model proposed in Ref. 1 is of the form 


T(t) = Ty exp (- pat ) (4-1) 


where 


= vortex circulation strength as a function of time t 
= initial vortex circulation strength at t=0 


= decay coefficient 
= turbulence level 


mo oO 
el 


= vortex pair semispan 


This expression was suggested as a possible model for vortex decay in the atmosphere in 
Ref. 5, with a decay coefficient value estimated to be 0.41 . 


Again, a least squares approach was used to curve fit the vortex circulation Strength data 
developed by the generalized algorithm. The vortex circulation strength profiles were fit to Eq. (4- 
1) by minimizing the logarithmic error 


E = )[G-(A-Bt)/? (4-2) 


over all applicable time increments in each run examined, where 


G = logarithm of I; 

A = logarithm of fy 

B = bd/s 

I; = vortex circulation strength generalized algorithm values 


Equation (4-2) returns the best fit values for I, and B_ for each run examined 


(comparisons were shown in Figures 3-3 to 3-28 above). Since the semispan s was known, the 
value of B could be used to infer the product bq for each run examined. 


The product bq may be plotted as shown in Figure 4-1. Here the F-15 and F-16 data are 
compared with the previous F-4 data. The plot is given for flight speed. The average value of bq 
is 0.74 m/sec, with a standard deviation of 0.48 m/sec, consistent with Program WIND data 
(Ref. 1), but 50% higher in mean value. 


The first five seconds of each run can be averaged for the U and V horizontal wind 
velocities at the two U-V-W anemometers. This data is shown in Figures 4-2 to 4-3. It may be 
used to estimate the turbulence level q_ by a least squares curve fit to a neutral logarithmic profile 


of the form 
U = U,gin( zz) (4-3) 
where 
U = _ horizontal surface velocity 
Ug = reference velocity 
= reference height 
o = Surface roughness 


A consistent turbulence level may then be obtained by using an equation found in Ref. 2 


Joe 2 : 
q* = 0.845 U, (4-4) 


The turbulence results are given in Figure 4-4. 
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Figure 4-2: Horizontal velocities for October 15, 1990, at Dugway: U_ velocity normal to 


tower grid (solid); V velocity parallel to tower grid (dashed); tower 0 (top); 
tower 9 (bottom). 
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Figure 4-3: Velocities for October 26, 1990, at Dugway. 
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Figure 4-4: Turbulence levels at Dugway. 


4-6 


5. CONCLUSIONS 


Dugway high speed aircraft tower anemometer data has been shown to be consistent with 
previous tower data. This suggests that a decay parameter of 


bq = 0.56 m/sec (5-1) 


could be used for lower flight speeds, and 
bq = 0.74 m/sec. (5-2) 


for higher flight speeds. 
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